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We describe experimental signatures of Majorana fermion edge states, which form at the interface 
between a superconductor and the surface of a topological insulator. If a lead couples to the 
Majorana fermions through electron tunneling, the Majorana fermions induce resonant Andreev 
reflections from the lead to the grounded superconductor. The linear tunneling conductance is 
{2e 2 /h) if there is an even (odd) number of vortices in the superconductor. Similar resonance occurs 
for tunneling into the zero mode in the vortex core. We also study the current and noise of a 
two-lead device. 



Introduction — Verifying the existence of Majorana 
fermions in condensed matter systems is an important 
topic in recent years, because of their potential appli- 
cation for quantum computations which are free from 
decoherence. [TJ [2] Quantum Hall states as well as super- 
conductors and superfluids with p x + ip y pairing symme- 
try are candidates which support Majorana fermions. [2] 
HJ |31 03 [51 [7J [HI [H] However, Majorana fermions in those 
systems are yet to be found. 

Recently, Fu and Kane [10 proposed that Majo- 
rana fermions can be created in the vortices of s-wave 
superconductors deposited on the surface of a three- 
dimensional topological insulator. [TTJ [T2J [T3J HU [T5J [T5] 
Moreover, chiral Majorana fermion edge states can be 
created at the interface between a superconductor and 
the area gapped by ferromagnetic materials |10j and 
several experiments with rather complex geometry have 
been proposed to study them. [TTl ITS] In this work, we 
propose experiments with relatively simple geometry to 
probe the chiral Majorana fermion edge states. 

More specifically, we study the tunneling current and 
noise from non-interacting Fermi leads to a grounded 
superconductor which possesses chiral Majorana edge 
states at its boundary. The experimental setup is shown 
in Fig. 1. An s-wave superconducting island is deposited 
on the surface of a topological insulator whose surface 
state is described by gapless Dirac fermions. [TTJ [12] 
The area outside the superconductor is gapped by fer- 
romagnetic materials. At the interface between the su- 
perconductor and the ferromagnetic material, there are 
gapless chiral Majorana fermion modes surrounding the 
superconductor. [10l H7J HHj One or two non-interacting 
Fermi leads are coupled to this chiral Majorana modes 
at points a and b with amplitudes t\ and t 2 respectively. 

We first consider the single-lead case by setting t 2 to 
zero in Fig. 1. We show that Majorana fermions in- 
duce resonant Andreev reflections from the lead to the 
grounded superconductor and result in highly nonlinear 
I-V curve which is a set of steps. At small voltage, the 
conductance from the lead to the superconductor is 0. 
However, the presence of a vortex in the superconductor 
changes the conductance to 2e 2 /h. For a two-lead de- 




FIG. 1: A superconducting island is deposited on the surface 
of a three-dimensional topological insulator. The area outside 
the superconductor is gapped by ferromagnetic materials. At 
the interface between the superconductor and the ferromag- 
netic material, there is a branch of chiral Majorana fermions 
denoted by 77. Two non-interacting leads are coupled to the 
Majorana fermions at point a and point b with amplitudes ti 
and t% respectively. The superconductor is grounded. 

vice with non-zero t 2 , crossed Andreev reflections may 
happen. In the small voltage regime, crossed Andreev 
reflections dominate over local Andreev reflections and 
the cross current-current correlations of the two leads 
are maximally positively-correlated. On the other hand, 
the presence of a vortex in the superconductor in the 
small voltage regime increases the current dramatically 
and the cross current-current correlations become maxi- 
mally negatively-correlated. 

Single-lead Device — This two-terminal device is shown 
in Fig. 1 by setting t 2 to zero. For simplicity, let us 
assume that there is only a single mode in lead 1. The 
Hamiltonian of lead 1 is: 

r+00 

H L1 = -iVf ^2 ^2 *PLa( X )dxlplea(x)dx, (1) 

e=R,L cr=1,l J ° 

where the tip of the lead is located at x — and Vf 
denotes the Fermi velocity. Hl\ contains both left and 
right moving fields but setting ipiL^ix) = ipiB.a(-x) f° r 
x > maps the Hamiltonian into one described by chiral 
fields only. 

/•+0O 

H L i = -ivf ^2 / yjl a (x)d x yj 1(J (x)dx. (2) 
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Let Hmo describe the chiral Majorana fermion mode sur- 
rounding the superconductor [TS] and Hti the coupling 
term between the lead and the chiral Majorana fermion 
mode with coupling strength t\. The total Hamiltonian 
is: 



H 1 = H 



LI 



H 



A/0 



H 



where 



and 



H M o = iv m Jq r]{x)d x ri{x)dx, 
Hti = ~i^h £ CT=T , r,(a)[^ lrT (0) + CC(0)]- 



(3) 
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In Eq|4j v m denotes the Fermi velocity of the Majorana 
mode. It is important to note that the Majorana mode 
on the surface of the topological insulator includes both 
spin components and has four instead of two components 
as in the case of spin polarized p-wave superconductor. 
Thus, the Majorana mode in Eq|4] couples to both spin 
up and spin down electrons. £ CT are complex numbers 
with \&\ = 1. 

We define two Fermi fields: ipi = ^(^tV'iT + 
Cl^il) and = t^tV^it ~ £i^u) such that the 
Majorana fermion couples to ipi only. Dropping ip[ 
we have: H[ = H' L1 + Hmo + Hti where H' L1 = 
~ iv f J-tt ip\(x)d x il>i(x)dx and H T \ = -itir)(a)[ipi+ip{]. 

In order to calculate the tunneling current from the 
lead to the grounded superconductor through its chiral 
Majorana fermion modes, we first calculate the scatter- 
ing matrix of the model described by H[ . Denoting the 
incoming and outgoing scattering states of the electrons 
and holes by Vifc(O-F) and ip\_ k (0^f) respectively, the 
scattering matrix can be written as: 



Vhfc(0+) 

v4- fe (o+) 



= s 



^ifc(O-) 



(5) 



where 



S 



1 

z 



zsin(0/2) -t?cos(0/2) 
-??cos(0/2) isin(0/2) 



(6) 

In Eq|6] Z = isin(0(fe,n)/2) + t%cos(9(k,n)/2) and 
9(k, n) = kL + 7r + nix is the phase a Majorana fermion 
with wave vector k acquires when it makes a complete 
circle around the superconducting island. L is the cir- 
cumference of the island and tt is Berry phase contribu- 
tion from the spin, n is the number of vortices in the 
superconductor. t% = t\/{2^/v m Vf) is dimensionless. 

From the structure of the scattering matrix, it is obvi- 
ous that when (9(fc, 0)/2 = kL + ir = 2mn with integer to, 
\s | 2 = 1. In other words, when an incoming electron 
has energy which matches an energy level of the quan- 
tized chiral Majorana modes, there is a resonant Andreev 
reflection. This means that an incident electron from the 
lead is converted into a backscattered hole with proba- 
bility of unity, independent of the coupling strength. We 





FIG. 2: a) Conventional resonant tunneling: two leads are 
coupled to a resonant level with coupling amplitudes ti and 
t2 respectively. Resonant tunneling with unity transmission 
probability can happen only if ti = t%. b) MIRAR: a single 
lead coupled to a Majorana level plays the role of both an 
electron lead and a hole lead. The coupling amplitudes of 
the leads to the Majorana level are ensured to be the same. 
The Majorana mode is attached to a superconductor which is 
grounded. 



call this Majorana fermion induced resonant Andreev re- 
flection (MIRAR). This is in sharp contrast to the usual 
normal metal-insulator-supcrconductor junction in which 
the Andreev reflection amplitude at fixed subgap energy 
always decreases with decreasing coupling strength. 

To acquire a physical picture of MIRAR, we note that 
for conventional resonant tunneling, unity transmission 
requires tuning the coupling strengths of the resonant 
level with the two leads to equal amplitude (Fig. 2a). 
However, when a lead is coupled to a Majorana fermion 
mode, the lead plays the role of both an electron lead and 
a hole lead. Because of its self-Hermitian property, the 
Majorana fermion is ensured to couple to the electron and 
hole leads with equal amplitude as can be seen in Eq|4] 
This results in resonant tunneling from the electron lead 
to the hole lead (Fig. 2b). 

The physical picture shows that MIRAR is a very gen- 
eral phenomenon which happens whenever a discrete Ma- 
jorana state is coupled to a Fermi lead. Another interest- 
ing example is the tunneling into the Majorana zero mode 
expected to exist in the vortex core. This experiment can 
be performed using a STM tip. The predicted linear con- 
ductance of 2e 2 /h should be a spectacular signature of 
the zero mode. When two vortices approach each other, 
the zero modes are coupled and split into occupied and 
unoccupied fermions and the resonance disappears once 
the energy splitting exceeds the voltage and the temper- 
ature. The mathematics is similar to tunneling into two 
localized Majorana states considered by Nilsson et aZ.,[T5] 
even though they focus on the regime of large splitting 
and did not consider the resonant tunneling into an iso- 
lated vortex. 

From the scattering matrix Eq(6j the tunneling current 
from the lead to the superconductor is: 



7 = 



2e 



eV 



T(E)dE 



2c 



eV 



ff(l + cos0) 



1 -cos0-Mf(l 



cos 9) 
(7) 

where T(E) denotes the Andreev reflection probability 
\s he \ 2 at energy E. Near cos0 = 1, T(E) can be cast 



dE. 
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FIG. 3: dl/dV verses eV with t\ = 0.1. eV is in units of 
■Khv m /L and dl/dV is in units of Solid (dashed) line 

represents the case with even (odd) number of vortices in the 
superconductor. 



into the resonance form: T(E) = 
7i = 2t\v m h/L and E t = (21 



7i 



where 



(B-B,) 2 +7 1 2; 

i = {M ■+■ n + 1) hv m Tt / L denote the 
quantized energy levels of the chiral Majorana fermion 
modes. Provided 71 is less than the level spacing, the 
differential conductance dl/dV verses eV peaks at 2e 2 /h 
whenever the electron energy is resonant with the Majo- 
rana mode energy Eg. As seen from Fig. 3 the resonance is 
shifted in voltage by half the level spacing when a vortex 
is added. For eV -C 7, the conductance jumps between 
2e 2 /h and near zero. We consider this a clear signature 
of the Majorana mode. 

Two-lead Device — Next we couple to the chiral Ma- 
jorana mode with one more lead, by setting t% to non- 
zero in Fig. 1. The new Hamiltonian becomes: H 2 = 
H[ + H' L2 + Ht2, where H' L2 and Ht2 are H' L1 and Hti 
with -0!, t\ and a replaced by ip 2 , ^2 and b respectively. 

The scattering matrix in the basics of 
(i/'ifc;'/'2fc,i/ , i_fc' can be written as: 



S = 



1 + A A 
A 1 + A 



(8) 



where 



Z' 



itjtj sin(§) - t{ cos(f) 1^*1*2 cos(|) 



1^*1*2 cos(|) iiltlsin(|)-ilcos(|) 

(9) 



In the above equation, Z' = sin(0/2) + (i?+ 

t 2 ) cos(0/2), ti — t i /(2 y /VfV m ). 9(k, n) = kL + n + rnr is 
the same as in the two terminal case. a((3) is the phase 
factor acquired by a Majorana mode propagating from 
point a(b) to point 6(a). 

Because of the special form of the scattering matrix, 
the average current from lead i to the grounded super- 
conductor /, 

be written in a compact form: 



and the current noise correlators P,j can 



Pij eliSij 



2f 
h 



2e 

reV 



eV 



(10) 



f [\A ij + (AA%\ 2 -\(AA%\^dE, 
Jo 



where the current noise correlators are defined as 



Pi, 



+ OC 



< [7,(0) -Ii][Ii(t) 



(12) 



The total current from the leads to the superconductor 
I = h + h is 



2e 



eV 



where 



T(E) = 1 



T(E)dE, 



(1-^)8^(0/2) 



(13) 



(1 



i 2 t 2 2 ) 2 sin 2 8/2 - 



{t 2 + t 2 ) 2 cos 2 (9/2) 



(14) 

For t\,t 2 < 1, the resonant Andreev reflection condition 
is E — Ei « 2hvf(t 2 + t 2 )/L, where Ei are the energy 
levels of the Majorana modes. As in the two-terminal de- 
vice, the I-V curve of the total current is highly non- linear 
and is a set of steps. In the small voltage regime with 
eV « 2hvf(tl + t\)/L, the conductance is ^~l~ 2 ^ 
when there arc no vortices in the superconductor. It is 
important to note that the linear conductance is not as 
in the two-terminal case. We argue below that this is the 
consequence of crossed Andreev reflections. If a vortex is 
created in the superconductor, the conductance becomes 
#- because of resonant Andreev reflection induced by 
the Ei — Majorana mode. 

Crossed Andreev reflection is a process which an in- 
coming electron from say, lead 1, is turned into an out- 
going hole in lead 2. As a result, one electron from each 
lead tunnels into the superconductor to form a Copper 
pair. This process is not allowed in the single-lead de- 
vice. We show below that measuring the shot noise of 
the tunneling currents can be used to reveal the mecha- 
nism of the tunneling processes, whether they are due to 
local or crossed Andreev reflections. In the following, we 
analyses the shot noise in the small voltage regime with 
■eV « 2hv f (i\ + i 2 2 )/L. 

In this regime and in the absence of vortices, we have: 



h=h = 1/2 = 



2e 2 V 



L x l 2 



h 1 



t\t\ 



(15) 



It is important to note that the individual contributions 
to the total current I from the two leads are equal and 
depend on the products of the coupling strengths t\t2 
only. This is a strong indication that the tunneling pro- 
cesses are dominated by crossed Andreev reflections. In 
order to verify this, we note that the noise correlators are 
given by: 



Pli — P22 — Pi 2 — P21 



2e 3 V 



72X2 
L l L 2 



(16) 



From Eq 



16 



we have Pa 



h (l + t 2 t 2 ) 2 
fgsli. For (ht 2 ) 2 « 1, 



(11) the Fano factor which is defined as Pu/eli approaches 1 
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indicating that crossed Andreev reflections dominate in 
this small voltage regime. 

Another point is that the cross current-current corre- 
lation function P 12 = (Pn + P22V2. As pointed out 
in Ref|T5} for any stochastic process the cross correlator 
has to satisfy the condition: |Pi 2 | < §(Pli + p22)- Thus, 
the currents from the leads are positively correlated in a 
maximal way, which is a consequence of crossed Andreev 
reflections and consistent with the finding of ReffTSIwhere 
two leads are coupled to two localized Majorana modes. 

In the small voltage regime with odd number of vor- 
tices, the average current from lead % to the superconduc- 



tor is / - *' 



and the total current is / = 2e . v , 



independent of the coupling strengths. In this regime, 



P11 — P22 — 



t\l\e 2e 2 V 



(17) 



A few comments of this result are in order. First, the 

noise power of each individual leads are non-zero and 

P - P - 

Pli = Israeli = P 22 = p^j|el 2 . To understand this 

result better, let us assume t± << t 2 . In this limit, 
Pli ~ eli. In other words, the Fano factor is 1. Phys- 
ically, it means that the tunneling events in lead 1 are 
dominated by crossed Andreev reflections, which give a 
Fano factor of 1 instead of 2 as in the case of tunneling 
into a conventional superconductor which is dominated 
by local Andreev reflections. 

On the other hand, the Fano factor at lead 2 ap- 
proaches zero. This is expected because the majority of 
the current from the two leads to the superconductor is 
carried by lead 2. When ti «t 2 almost all the incoming 
electrons from lead 2 are locally Andreev reflected. This 
is analogous to the situation of a ballistic tunneling junc- 
tion with tunneling probability T where the Fano factor 
is suppressed by a factor of (1 — T). For T « 1, the Fano 
factor approaches zero. 

One more interesting point is that the total noise power 

Pr = 1^=1,2 A? is zero and p i2 = -|(Pn + p ii)- In 
other words, the two currents are negatively correlated in 
a maximal way. It can be argued that this is a result of 
crossed Andreev reflection and MIRAR. 

Conclusion and Discussion — We show that Majorana 
modes which couple to Fermi leads induce resonant An- 
dreev reflections from the leads to the superconductor. 
At small voltage limit, the conductance from a single 
lead to the superconductor is (2e 2 /h) if there is an even 
(odd) number of vortices in the superconductor. Attach- 
ing one more lead to the Majorana fermion edge states 
introduces crossed Andreev reflection. The currents from 
the two leads can be maximally positively-correlated or 
maximally negatively-correlated, depending on the par- 
ity of the number of vortices. 

Our discussion can easily be extended to finite tem- 
perature T and it is clear that the effects we discussed 
will be smeared out if the temperature is larger than 



the level spacing of the Majorana fermions, T > Tiv m /L. 
This puts a rather stringent, but not impossible, require- 
ment on the size of the dot and position of the Fermi 
level relative to the Dirac node in order to minimize L 
and maximize t>m.|18] On the other hand, resonant tun- 
neling into the Majorana zero mode in the vortex core 
is subject to the less stringent condition T < AE where 
AE ~ A 2 /Ep is the level spacing of states in the core. 
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